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ABSTRACT 

We examine recent evidence from the luminosity-redshift relation of Type la Supernovae (SNe la) for 
the ~ 3 a detection of a "H ubble bubble" - a departure of the local value of the Hubble constant from 
its globally averaged value (|Jha et al.l I2007T ) . By comparing the MLCS2k2 fits used in that study to 
the results from other light-curve fitters applied to the same data, we demonstrate that this is related 
to the interpretation of SN color excesses (after correction for a light-curve shape-color relation) and 
the presence of a color gradient across the local sample. If the slope of the linear relation {(3) between 
SN color excess and luminosity is fit empirically, then the bubble disappears. If, on the other hand, the 
color excess arises purely from Milky- Way like dust, then SN data clearly favors a Hubble bubble. We 
find that SN data give /3 ~ 2, instead of the (3 ~ 4 one would expect from purely Milky Way- like dust. 
This suggests that either SN intrinsic colors are more complicated than can be described with a single 
light-curve shape parameter, or that dust around SN is unusual. Disentangling these possibilities is 
both a challenge and an opportunity for large-survey SN la cosmology. 
Subject headings: cosmology: observations — supernovae: general 



1. INTRODUCTION 

In an analysis of th e luminosi t y dist ances of Type la 
supernovae (SNe la), Uha et all (12001 (hereafter J07) 
presented evidence for an offset between the Hubble con- 
stant measured from SNe la with 2500 < cz < 7400 km 
s _1 and from those with 7400 < cz < 45000 km s" 1 . 
Specifically, the more distant SNe are slightly fainter 
than one would expect, implying a high local value of 
Ho. One natural explanation is a "Hubble bubble" - 
a local monopole in the peculiar velocity field, perhaps 
caused by a local void in the mass density. The analysis 
of J07 was carried out using the most recent manifesta- 
tion of the MLCS light-curve fitter, MLCS2k 2, and found 
SH/H = 6.5% ± 1.8%. The Am 15 method dPrieto et all 
gives similar results. In c ontrast, galaxy-cluster 
distances give SH/H ~ 1.5%±2% (|Giovanelli et al.lll999fc 
iHudson et all I2004D . which is marginally inconsistent 
with the SN result (1.9&-). E a rlier s tudies using SNe la 
can b e found in iKim et ail {T997) and IZehavi et al.1 
(|1998ft . A Hubble bubble could have serious implications 
for precision SN la cosmology programs. 

In this Letter we test the evidence for the Hubble bub- 
ble using three ot her light-c urve analysis pac kages: SALT 
(|Guv et alJl2005h . SALT2 (|Guv et al.ll200l . and an un- 
published package dev eloped for the 3rd year Supernova 
Legacy Survey fSNLS. lAstier et al.ll20M ) data (SiFTO). 
An analysis of the same data set with these tools does 
not support a bubble if SN data are used to derive the 
relationship between SN color excess and peak luminos- 
ity, but does if this relationship is required to be that of 
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Milky Way- like dust. Therefore the question of the Hub- 
ble bubble is one of the nature of SN co l ors. A similar 
result was found independently bv I Wand (|2007l ). 

2. LIGHT-CURVE FITTERS AND THE MEANING 
OF SN COLORS 

A considerable amount of effort has been devoted to 
the question of how to best fit SN la light curves. The 
goal is to measure the relative luminosity distance of 
different SNe, usually after correcting for the width- 
luminosity and color-luminosity relationships. These ap- 
proaches differ considerably in implementation and as- 
sumptions, but generally produce very similar results 
when used to e stimate the cosmological density parame- 
ters (O m ,f2 A ) (|Wood-Vasev et al.ll2007T) . 

In this analysis we consider four fitting packages: 
MLCS2k2, as well as three packages developed for use 
with SNLS: SiFTO, SALT, and SALT2. The latter three 
are more similar to each other than to MLCS2k2. The 
details of SiFTO will be presented elsewhere (Sullivan 
et al., in preparation); it is broadly similar to SALT ex- 
cept that the color excess relation is not imposed during 
the light-curve fit, but rather when the results of the fit 
are converted into a distance estimate. 

The critical issue for the question of the Hubble bubble 
relates to how the different fitters handle SN color ex- 
cesses. All of the models have two parameters: one that 
describes the light-curve shape, and one the color, both 
of which affect the peak luminosity. For MLCS2k2 these 
are A and Ay, for SALT they are s and c. These have 
different technical meanings, but essentially describe the 
same effects. MLCS2k2 breaks the color relation into two 
pieces: a time- independent piece 6 , parameterized by Ay, 
and a time-dependent piece determined by the shape of 
the light curve, A. The t ime independent pi ece is as- 
sumed to follow the CCM (|Cardelli et al.lll989D dust ex- 

6 Technically, the effective value of Ry depends on the SED, so 
the effects of this term are weakly time-dependent. 
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tinction law. In effect, Ay parameterizes the color of the 
SN after correction for a shape-color relation. SALT fol- 
lows a similar prescription, with some critical differences. 
The wavelength behavior of the time-independent piece 
(parameterized by c) is derived as part of the SALT train - 
ing process; the results are shown in iGuv et al.1 {2007). 
This relation is reasonably similar to the CCM law, and 
is arbitrarily normalized so that c is the B — V color 
at peak luminosity (plus a small constant) for any s; the 
same is not true of U — B and V — R, which are functions 
of s. In the following we refer to the time-independent 
piece as the color excess relation. 

The real difference is how the color excess relation is 
used to correct the peak luminosities. MLCS2k2 makes 
the simplest assumption - that it is due to dust. The 
wavelength behavior of the CCM law is parameterized 
by the ratio of selective to total extinction (Rv), and 
this is also used to convert the color excess into the mag- 
nitude correction. The other fitters considered here take 
a different approach. One of the steps in the process of 
turning their results into a distance estimate is to em- 
pirically measure the slope of the relationship between 
the color c and the luminosity correction ((3), as well 
as between the light-curve shape and luminosity (a) by 
minimizing the residuals with respect to the Hubble line. 
Here the model for the predicted peak B magnitude m b 
is m B = Mb + 5 log 10 d^ — a (s — l)+/3c, where M B is the 
absolute peak magnitude, dj, is the luminosity distance, 
and a and (3 are determined from SN data; the slopes 
can be cleanly separated from the intercept without any 
difficulties. In this formulation, the shape-color relation 
is absorbed into the a coefficient, and therefore the (3 
correction only applies to the color excess. In this analy- 
sis we generally work with B magnitudes, and therefore 
if the color excess relation is purely due to dust then we 
expect (3 — Rb = Rv + 1 — 4. Interestingly, the current 
constraint on (3 from combined SNLS and low-z data is 
j3 ~ 2 ± 0.2, which differ s from 4 by > 10<r. A similar 
value was found by EES (EMI). 

What does j3 ^ 4 mean? If (3 is interpreted as aris- 
ing from dust, then this would require Rb = 2, which is 
extreme. On the other hand, we know that dust must 
be present at some level. One possible explanation is 
that the various color parameters are actually measuring 
some combination of two effects: dust, as well as some 
additional intrinsic color variations that are not related 
to A or s. The (3 for this additional piece is not known, 
but if it is less than the value for dust, then the effec- 
tive combined value will be driven down. This additional 
color parameter could be interpreted as additional scat- 
ter in the measured colors of SN la; most of the fitters 
studied here (except SALT) do include terms for such 
scatter, but assume that it has no effect on the overall 
luminosity ((3 = 0). A worry, then, is that the relative 
balance of dust and the additional intrinsic color might 
vary between SNe in different environments. Alterna- 
tively, it is possible that the dust along the line of sight 
to some SNe la is different than standar d dust, perh aps 
due to scattering effects local to the SN (Wang 2005). 

3. TESTING THE HUBBLE BUBBLE WITH 
OTHER FITTERS 

The data sample used in the analysis of J07 incorpo- 
rates a number of SN with sparsely sampled light curves. 
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Fig. 1. — SALT2 c parameter compared with the equivalent 
MLCS2k2 value, A v /R v for our 61 SNe. The A v prior has been 
removed from the MLCS2k2 fits to allow easier comparison of 
the underlying fits. Good agreement is also found in the times 
of maximum, peak magnitudes, and light-curve shape parameters 
(although for the latter the relations are non-linear). 

We have placed the following additional requirements on 
our sample: First we require at least one rest frame B 
observation within ±7 days of the estimated B peak, in 
order to ensure that tub is well measured, and at least 
one U or V observation within -7 to +10 days for the 
color. We exclude all SN with s < 0.7 or s > 1.3, since 
these are outside the model bounds for most of the fit- 
ters. Finally, we require that E ( B — V) MW < . 4 mag 
as calculated from the dust maps of lSchlegel et al.l (|1998f ) 
and c < 0.6 mag, both to avoid regions of anomalous 
dust. This reduces the sample of 95 SN used by J07 to 
61 objects. Since only MLCS2k2 is trained to fit / data, 
we remove this from all fits. In order to check that this 
reduction in the sample does not affect the results, we 
then refit all of these SN using the most recent version 
of MLCS2k2. We find good agreement with J07, with 
SH/H = 6.0% ± 1.9% for cz yoid = 7400 km s" 1 and 
assuming a flat fi m = 0.3 ACDM universe. 

We then fit the same photometry using the other light- 
curve packages 7 . The different packages show impres- 
sive agreement in their basic derived light-curve param- 
eters, as demonstrated in figure [TJ In order to measure 
SH/Hwith SALT/SALT2/SiFTO, we must determine (3. 
This is calculated from the SN data, but incorporating 
high-redshift SNLS data as well. SALT, SALT2, and 
SiFTO give (3 = 1.82, 1.75 and 2.31, respectively, with 
errors of about ±0.16. The different values are not par- 
ticularly important for the current analysis; all three 
can be fixed to the mean value without qualitatively 
changing the results. After applying the color and light- 
curve shape corrections, we obtain SH/H = 0.9%±2.0%, 
-1.2% ± 2.1%, and 0.4% ± 2.1% - in other words, they 
do not support a Hubble bubble if (3 is fit to the data. 
SALT and MLCS2k2 are compared in figure [2] for var- 
ious values of cz VQ ^. The exact values are somewhat 
sensitive to the cuts applied to the sample. However, it 
requires considerable hand tuning to find regions of pa- 
rameter space where any of the fitters besides MLCS2k2 
detect a Hubble bubble at greater than 1.4 a, and in all 
cases MLCS2k2 gives much greater significance. 

One might be tempted to interpret this as an er- 

7 The raw light-curve parameters are available from 
http : //qold. astro .utoronto . ca/conley/bubble 
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Fig. 2. — Hubble bubble as a function of the velocity of the 
step in the Hubble constant, cz VQ ^. Left: Results for SALT (with 
P = 1.82). Right: Results for MLCS2k2. The grey band is the 
error in SH/H. The constant value for 16000 < cz < 19000 km 
s _1 simply reflects the lack of any SN in this range. 




Fig. 3. — SH/H vs. the value of /3 used to convert the measured 
color parameter into a luminosity correction for SiFTO. The rela- 
tion is almost, but not quite, linear. Fits to low-z and SNLS data 
give /3 = 2.35 ± 0.16. 



ror in one or more of the packages. However, the ac- 
tual cause is more interesting: if one sets j3 = 4.1 for 
SALT/SALT2/SiFTO, they do find evidence for a Hub- 
ble bubble at > 2.5<j, as shown in figure [31 The color 
model for the most recent version of Amis, which also 
favors the bubble, is similar to that of MLCS2k2. The 
question of the Hubble bubble therefore boils down to the 
appropriate value of (3 - is it the ~ 4 of Milky Way-like 
dust, or is something more complicated going on? 

The value of (3 is relevant because the nearby portion 
of the low-z SN sample is redder than the distant por- 
tion (figure |H). This is probably due to Malmquist bias 
(|Mahnquist et al] [T936) or other selection effects, since 
redder SNe are fainter and harder to detect. The effects 
of the color correction are to make the blue SNe dimmer 
and the red SNe brighter (relatively). Malmquist bias 
has little effect as long as the appropriate value of (3 is 
applied across the whole sample; however, if the right 
value is j3 = 2 and instead 4 is used, then the distant 
portion will be made too faint, and the nearby portion 
too bright, which is exactly the effect observed in J07. 
MLCS2k2 usually includes a prior on Ay, and the default 
version does not take into account the redshift-dependent 
effects of Malmquist bias. Adjusting the prior to reflect 
this increases the value of SH/H. 
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Fig. 4. — Colors of the low-z supernova sample using SALT. The 
distant portion of the sample is bluer than the nearby portion, 
which is probably caused by selection effects. 

4. THE VALUE OF (3 

If, in fact, (3 = Rb (i-e., it is purely dust), then we 
should note that a large range of Rb values have been 
observed in different environments. However, a sam- 
ple mea n value of Rb = 2 would be quite surprising 
praind l2003h . even with selection effects. Most pre- 
vious studies of large samples of SNe la have found 
R B ^ 3.5 ± 0.3 (Idella Valle fc Panagialll992t IRiess et all 
[19911 iPhillms etal.1 119991 : lAltavilla et alJ [2004]). This 
value weakens the evidence for the Hubble bubble slightly 
to ~ 2 a. 

Our empirical values for (3 were calculated in a model 
that assumed a smooth local Hubble flow. It is possible 
that this could be artificially suppressing 8H/H, so we 
checked this by refitting (3 using only high-z SNLS SN, 
only the low-z SN below cz yQ ^ and simultaneously with 
our fits to SH/H. These give essentially the same result 
(j3 ~ 2), albeit with larger errors than the full sample. 
Therefore, (3 is robust against the presence of a Hubble 
bubble. Monte Carlo studies indicate that the bias in (3 
{bp) due to errors in the measurement uncertainties and 
covariances is \bp\ < 0.02 for fairly extreme cases. 

We can also analyze the results of the MLCS2k2 fits 
using the (3 framework. This differs from what is meant 
in J07 by fitting Ry , which enforces a certain relation be- 
tween the wavelength dependence and scaling of the color 
excess relation. We first remove the extinction correction 
from the MLCS2k2 distance estimates, then convert Ay 
into E (B — V) using Ry, and finally use this to fit for 
the value of (3 by minimizing the residuals with respect 
for the Hubble line via a (3E (B — V) term. We work in 
B, and find = 2.7 ± 0.3, which again should be com- 
pared with the expected value of 4.1. Restricting the fit 
to only SNe below cz v0 ^ does not change the results. 
The fits are shown in figure 

Another technique for estimating the amount of ex- 
tinction is to use late time (~ 45 days after peak) 
color measurements. The idea is that at late times all 
SNe la have a simple relationship between intrinsic color 
and epoch, and this can be used to measure extinction 
([Phillips et al.lll999L hereafter P99). The evidence for 
this is based on a handful of SN for which there is inde- 
pendent evidence for low extinction, and the distribution 
of the measured late-time colors (J07, figure 6). A ver- 
sion of this is used in the training process of MLCS2k2. 

This suggests one more test of (3. We take the subsam- 
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Fig. 5.— Fits to for SiFTO (top left), SALT (top right), 
and MLCS2k2 (bottom left). The residuals are compared to the 
best fitting Hubble line but without correction for the various color 
parameters, which are shown along the abscissas. Bottom right: 
Results for SiFTO using the late time colors as described in the 
text. The solid lines are the best fit to /3 as given in the text. 
The dashed lines show /3 = 4.1, which is the expected value if the 
color excess relation is caused by Milky Way-like dust. The black 
circles are for SNe with cz < 7400 km s —1 , and the blue squares 
for cz > 7400 km s _1 . The grouping around Ay/Ry = for 
MLCS2k2 is the result of the Ay prior. 

pie of our 61 SNe that have late-time color measurements 
(from J07) and use these values as our color estimate 
to fit for (3. Carrying out this analysis out for all four 
fitters, we find (3 ~ 2.3 ± 0.3 for 37 SNe. In order to 
eliminate any "cross-talk" between the color parameters 
and the peak magnitudes, we also tested this procedure 
using only B band data to fit the peak magnitude and 
light-curve shape and obtained the same results. Using 
the late-time colors provided by P99 gives /3 ~ 1.5 for 
28 SNe. These results again differ considerably from the 
expectation of 4.1. 

5. CONCLUSIONS AND DISCUSSION 

We have demonstrated that the SN la evidence for a 
Hubble bubble is related to how SN colors are modeled. 
All of the approaches considered agree that the wave- 
length dependence of the color excess relation is similar 
to that of dust, but disagree on whether or not the re- 
lation between the measured color excess and the peak 
luminosity is also dust-like. Our fits give a value of (3 ~ 2, 
which if interpreted via the CCM dust law, requires the 
extreme value Ry ~ 1. Therefore, either a more com- 
plicated model of intrinsic SN colors is required, which 
goes beyond a single light-curve shape-color relation, or 
dust in the host galaxies of SNe la is quite atypical of 
Milky Way dust. If the former, then the late-time colors 
of SNe la vary from SN to SN, and therefore do not pro- 
vide a simple measure of extinction. If one does favor the 
single-parameter model with Galactic dust, then the evi- 
dence for the Hubble bubble from SNe la is fairly strong. 



These results depend on how accurately the light-curve 
shape-intrinsic color relationship has been modeled, so 
it is reassuring that the four fitters are in approximate 
agreement on the value of j3. 

Requiring an additional intrinsic color relation beyond 
the shape-color relation raises some challenges for SN 
research. Unless this relation can be disentangled from 
dust, we must consider the possibility that the balance 
of the two effects will change with environment and red- 
shift, perhaps even within the low-z sample, and affect 
precision SN cosmology. MLCS2k2 might require a more 
sophisticated prior that takes into account the relation 
between extinction, light-curve shape, intrinsic color, and 
would also have to allow for any effects of this intrinsic 
color on SN luminosity. SALT/SALT2/SiFTO might re- 
quire different values of (3 in different environments or for 
different color thresholds. Using an inappropriate value 
of (3 will mostly affect the most distant SN in any sur- 
vey, where Malmquist bias is important. The potential 
systematic for a given SN survey can be evaluated by 
multiplying the uncertainty in f3 by the change in color 
excess across the sample. If the two-component color 
model is correct, then constraining it will not be a triv- 
ial task, but it does hold out the possibility of making 
SNe la even better standard candles. 

How can we resolve this issue? A deeper nearby su- 
pernova sample that has a similar color distribution at 



all distances out to 



25000 km s 1 would provide 



a good test of the Hubble bubble, since an incorrect 
value of (3 would no longer introduce such an effect. Note 
that SN cosmology an alyses that re s trict t h emselves to 
fmin > 0-015 such as lAstier et al.l (|2006[ ); iRiess et all 
(2007); IWood-Vasev et al.l <|200"7f T~i are not strongly af- 
fected by the existence of the bubble. Determining if 
a more complicated color model is necessary requires 
a different approach. A larger sample of SNe in low- 
extinction environments (like elliptical galaxies), or at 
least in a narrow color range, could be used to search for 
a non-dust-like color relation. A wider baseline of color 
measurements could also help this problem; it seems un- 
likely that the wavelength dependence of the color excess 
relationship will continue to look like dust at all wave- 
lengths unless it rea lly is dust. The S ALT/SALT2 color 
excess relationship (jGuy et al.l I2007L figure 3) displays 
tantalizing hints of departures from the CCM law. If 
these can be conclusively demonstrated, it would at least 
prove that there is more going on than Milky Way-like 
dust, even if it might not elucidate the underlying mech- 
anism. 



We thank Reynald Pain and Kathy Perrett for useful 
discussions. 
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